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Chapta*! 

INTRODUCTION 


In  recent  years  there  has  been  a  new  interest  in  millimeter  (mm)  wave 
technology.  Vacuum  tube  sources  that  cover  the  entire  mm  wave  spectrum  exist  but 
there  is  a  demand  for  solid  state  devices  capable  of  operation  at  mm  and  sub  mm  wave 
frequencies.  There  are  solid  state  sources  such  as  IMP  ATT  and  GUNN  diodes  that  can 
operate  in  the  lower  mm  wave  region  but  their  efficiency  and  performance  degrade  at 
higher  frequencies. 

Quantum  well  devices  have  recently  been  shown  to  have  a  potential  to  generate 
mm  and  sub  mm  wave  frequencies  [1].  This  is  because  the  intrinsic  negative  resistance 
region  for  the  quantum  well  device  extends  from  dc  up  to  several  hundred  GHz.  One 
drawback  of  the  quantum  well  device  is  its  low  output  power.  A  new  device,  the 
quantum  well  injection  and  transit  time  (QWTTT)  diode  was  proposed.  This  device 
exploits  transit  time  phenomena  to  improve  the  output  power  of  the  conventional 
quantum  well  devices.  Quantum  well  oscillators  are  also  very  attractive  because  they 
are  intrinsically  very  low  noise.  In  order  to  generate  higher  output  power  in  the  mm 
wave  region,  different  power  combining  techniques  can  be  used. 

Ordinarily  a  real  physical  system  is  a  nonlinear  function  of  its  operating 
parameters;  therefore,  if  a  complete  and  accurate  analysis  of  the  system  is  desirable  the 
nonlinear  equations  representing  the  physical  problem  should  be  solved.  Since  the  task 
of  solving  the  nonlinear  equations  is  often  very  difficult,  it  is  necessary  to  simplify  a 
nonlinear  problem  through  some  assumptions.  After  the  design  work  is  done  the 
system  can  be  built.  If  it  does  not  have  all  the  desirable  characteristics,  one  can  achieve 
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the  final  design  by  trial  and  error.  Depending  on  the  type  of  system,  the  designer  can 
impose  some  restrictions  on  the  operating  condition  of  that  system  so  then  its  operation 
in  a  linear  region  is  ensured.  In  this  case  small  signal  analysis  is  satisfactory. 

An  oscillator  relies  on  device  nonlinearity  for  proper  operation  and  is  inherently 
a  nonlinear  system.  Due  to  nonlinear  operation  often  higher  order  harmonics  are 
generated  that  effect  the  frequency  of  oscillation,  the  power  generated  and  the  efficiency 
of  the  oscillator.  It  is  only  by  nonlinear  analysis  that  one  can  determine  the  exact 
frequency  of  operation,  the  power  and  efficiency  of  the  oscillator,  the  power  of 
undesirable  higher  harmonics  relative  to  the  power  generated  at  the  fundamental 
frequency,  and  the  limitations  of  a  microwave  oscillator.  Another  important  fact  that 
makes  nonlinear  analysis  necessary  is  the  design  of  monolithic  microwave  integrated 
circuits  (MMIC).  The  opportunity  of  trimming  a  circuit  after  design  is  not  available  in 
the  case  of  MMICs  and  redesigning  and  reproducing  them  is  very  time  consuming  and 
costly.  Therefore,  it  is  very  desirable  to  be  able  to  design  and  construct  an  MMIC 
circuit  which  meets  all  the  specifications  without  going  through  any  trial  and  error.  The 
nonlinear  analysis  of  a  circuit  helps  the  design  engineer  to  achieve  this  goal. 

A  complete  analysis  of  a  nonlinear  circuit  requires  a  nonlinear  device  model  and 
analytic  means  to  extract  the  effect  of  device-circuit  interactions  from  the  model.  There 
exist  different  methods  to  solve  nonlinear  problems.  In  general  one  can  write  a  set  of 
nonlinear  differential  equations  to  represent  a  dynamical  system.  In  principle,  the 
analysis  can  be  performed  by  integrating  the  system's  differential  equations.  If  it  is 
desired  to  find  the  steady  state  response  of  a  circuit  as  it  is  in  the  case  osc 'iators,  it  is 
necessary  to  carry  out  the  integration  dll  the  transient  responses  die  out.  This  approach 
is  prohibitively  expensive  whenever  the  transients  are  governed  by  time  constants  that 
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are  very  much  larger  than  the  period  of  the  driving  force  of  an  amplifier  or  the 
oscillation  of  an  oscillator.  In  this  case,  continuous  integration  may  have  to  be  carried 
out  for  hundreds  of  full  periods  before  the  transients  become  small  enough  to  ignore. 
A  better  method  to  analyze  this  type  of  problems  is  the  harmonic  balance  technique  [2]. 
In  the  harm.onic  balance  method,  each  state  variable  is  represented  by  a  Fourier  series 
that  satisfies  the  requirement  of  periodicity.  An  optimization  algorithm  is  then  used  to 
adjust  the  coefficients  of  the  Fourier  series  such  that  the  system  equations  are  satisfied 
with  least  square  error. 

In  this  work  a  small  signal  and  large  signal  analysis  of  a  single  QWTTT  diode 
and  also  a  small  signal  analysis  of  a  power  combining  structure  are  performed.  At  the 
end  some  experimental  work  is  done  and  experimental  and  theoretical  data  are 
compared. 


Chapta‘2 

SMALL  SIGNAL  ANALYSIS 


In  order  to  realize  a  planar  QWITT  diode  oscillator  one  can  insen  a  QWITT 
diode  via  an  inductive  post  inside  a  parallel-plate  waveguide.  This  type  of 
configuration  is  convenient  for  design  of  monolithic  microwave  integrated  circuits 
(MMIC)  in  the  mm  wave  region. 

To  obtain  higher  output  power  at  mm  wave  frequencies  a  TEM  parallel-plate 
waveguide  can  be  loaded  periodically  with  QWITT  diodes.  This  simple  design 
includes  a  resonant  circuit  which  when  properly  designed  precludes  the  need  for  any 
external  resonant  circuitry. 

2,1  -  Analysis  of  The  Stnicture 

The  parallel-plate  waveguide  loaded  periodically  with  QWTTT  diodes  is  shown 
in  Fig.(2-1).  The  waveguide  is  open  on  one  end  and  terminated  by  a  resistive  load  at 
the  other  end.  Assuming  that  we  are  given  the  small  signal  impedance  of  the  QWETT 
diode  at  a  particular  frequency  (fr)  our  objective  is  to  design  a  periodic  QWITT  diode 
oscillator  to  oscillate  at  the  frequency  f,.  The  parameters  that  should  be  determined  are 
the  width  (w),  the  height  (b)  of  the  waveguide,  the  distance  between  the  load  and  the 
metallic  post  (din),  the  length  of  the  open  part  of  the  waveguide  (dout).  the  periodicity  of 
the  structure  (d)  (in  the  case  of  having  more  than  one  diode  mounted  inside  the 
waveguide)  and  the  impedance  of  the  load  (Zi). 
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Metal  Posts 


— 

IBI 

mm 
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QWnr  diode  or  Resonant 
tunneling  diode 

Fig.(2-1):  schematic  representation  showing  top  and  side  views  of  parallel 
plate  waveguide  periodically  loaded  by  QWTTT  diodes 


2.1.1  •  Calculation  of  The  Post  Impedance: 

The  parallel  plate  waveguide  loaded  periodically  with  the  QWITT  diodes  is 
modeled  by  an  equivalent  transmission  line  loaded  periodically  with  the  impedance  Zp. 
The  impedance  Zp  is  the  sum  of  the  small  signal  impedance  of  the  QWTTT  diode  and 
the  metal  post  reactance.  Here  we  treat  the  QWTTT  diode  as  a  lumped  element,  and  the 
electric  field  component  and  the  current  density  are  assumed  to  be  oriented  in  the  Y 
direction  (Fig.2-l(b)) 

The  small  signal  impedance  of  the  diode  and  the  metal-post  reactance  need  to  be 
calculated.  Since  the  thickne  ss  of  the  QWTTT  diode  is  much  smaller  than  the  length  of 
the  suppprting  metal  post,  we  can  not  neglect  the  effect  of  the  reactance  of  the  post. 


The  mode-matching  technique  can  be  used  to  calculate  the  first-order  inductance  of  the 
post  [3-4].  The  cross  sectional  view  of  the  waveguide  modeled  is  shown  in  Fig.(2-2). 
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Fig.(2-2);  Cross  sectional  view  of  the  analytical  model  for  the 
calculation  of  the  metal  post  reactance 


The  impedance  of  the  post  jX^,  can  be  obtained  by  the  following  formula. 


„  k  ,-a  .  ,  Tcr ,  .  ~  .  2m7iD  ,  ,  1  a  ,, 

Xo=j(-;^log(4sin(5j)sin(— ))+  Z  (  1  +  cos— 5— )  ( — -  — )]  (2-1) 

Ym 


Consider  that  X^,  is  the  normalized  reactance  of  the  post  with  respect  to  the 

characteristic  impedance  of  the  parallel-plate  waveguide.  The  above  formula  is  valid  if 

the  diameter  of  the  post  is  much  smaller  than  the  width  of  the  parallel  plate  waveguide. 
If  the  post  is  located  at  the  center  of  the  waveguide  (D  =  j),  then: 
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V  k,-a,  2nT  “ 1  a  ,, 

Xo  =  j{— log(-— )+  I2( - ^)} 

^  7C  *  m=2-4..  Y  mTC 

m 


where 


K=^f^K, 


The  total  post  impedance,  Zp,  is  the  characteristic  impedance  of  the  QWITT  and  the 
post  reactance  is  given  by  : 


2p  —  Z(jjQ(jg  +  JZqXq 


where  Z^,  is  the  characteristic  impedance  of  the  parallel-plate  waveguide.  Since  the 
parallel-plate  waveguide  is  modeled  as  a  TEM  transmission  line  then  referring  to 
Fig.(2-3)  we  can  write: 

V  =  E  b  (2-7) 


I  =  Ha 


I  “Ha~^I 


(2-10) 
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and 


Tio  =  120  K 


(2-11) 


I 

^  ^  z 

a 

Fig.(2-3):  TEM  model  of  a  parallel  plate  waveguide 

where  is  the  relative  dielectric  constant  of  the  medium  inside  the  parallel  plate 
waveguide.  The  above  formula  is  valid  when  the  width  of  the  parallel-plate  waveguide 
is  much  larger  than  its  height  (a»b). 

2.12  -  Small  Signal  Modeling  of  The  QWITT 

The  physical  structure  of  the  QWITT  diode  consists  of  a  single  GaAs  quantum 
well  sandwiched  between  two  thin  layers  of  Al^Caj.^As,  together  with  a  drift  region  of 
undoped  GaAs  [5].  This  structure  is  then  placed  between  two  n"*"  GaAs  regions  to 
form  contacts  (Fig.  2-4) 
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Fig.(2-4);  Phsical  structure  of  the  QWTTT  diode 


The  small  signal  equivalent  circuit  of  the  QWITT  is  shown  in  Fig.(2-5). 


Qrell 
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L>V\/V^ 


LJYYY) — VW— ' 
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prstc 
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Fig.(2-5);QWnTs  equivalent  circuit 


It  is  assumed  that  the  length  of  the  quantum  well  region  "1"  is  much  smaller  than 
the  depleted  spacer  layers  W.  The  specific  impedance  of  the  quantum  well  injection  at 
an  angular  frequency  co  is  then  given  by 

1 

a+jcoe' 


(2-12) 
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where  e'  is  the  effective  dielectric  constant  of  the  injection  region  and  a  is  the 
normalized  injection  conductance.  The  total  small  signal  specific  impedance  (excluding 
parasitic  elements)  of  the  QWTTT  diode  is  the  sum  of  this  specific  impedance  and  that 
of  the  depleted  spacer  region  Z^. 


Z  -  Zqw  +  Z 


n 


(2-13) 


Based  on  the  frequency-independent  injection  conductance  and  the  constant 
saturation  velocity  the  drift  region  electric  field  can  be  integrated  to  obtain  the  specific 
impedance  at  an  angular  frequency  O)  for  the  QWTTT  device  as 


W  j-^  g  l-exp(-jed)j 
jcoe  <T+jcoe  jO^ 


(2-14) 


where  £  is  the  dielectric  constant  of  the  drift  region  and  0^  is  the  drift  angle  given  by: 


(2-15) 


The  real  part  of  Eq.(2-14)  gives  the  negative  resistance  obtained  from  the  depleted 
spacer  region  as: 


(2-16) 


The  specific  resistance  of  the  quantum  well  injection  region  is  given  by  the  real  part  of 
Eq.(2-12) 


(2-17) 


The  total  specific  resistance  of  a  QWITT  diode  (excluding  parasitic  resistance  in  the 
device)  is  given  by  the  sum  of  Eqs  (2-16)  and  (2-17). 

2.U  -  Calculation  Of  Complete  Admittance 

Based  on  the  last  two  sections  we  can  determine  the  value  of  the  loaded 
impedance  Zp  which  is  the  sum  of  the  QWITT  diode  small  signal  impedance  and  the 
reactance  of  the  post.  The  periodic  structure  is  shown  in  Fig. (2-6). 


Fig.(2-6):  Transmission  line  model  of  the  finite  length  periodic  QV,  ITT 
oscillator 
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The  characteristics  of  this  structure  can  be  found  using  the  theory  of  periodic 
structures;  however,  a  different  approach  can  be  used.  In  practice  it  is  not  necessary  to 
have  many  diodes  and  it  is  desirable  to  keep  the  total  line  length  short  to  ensure 
operation  at  only  one  frequency.  That  is,  the  separation  of  the  resonant  frequencies 
becomes  large  as  the  length  decreases.  As  shown  in  Fig. (2-6),  simple  admittance 
transforms  can  be  used  to  analyze  the  periodic  soiicture. 

Y'l  +  tanh('ydout)  .  o 

Yj  _ - ; - +  Yp 

1+  YLtanh(ydout) 

y',  +  tanh(Yd)  . 

Y2  - - : - +  Yp 

1+  YjtanhCTd) 


Y'n  +  tanh(Ydin) 

1+  Y'„tanh(Ylin) 


Where  Y  is  the  impedance  normalized  to  the  charactristic  impedance  of  the  waveguide. 

If  we  cut  the  circuit  in  two  sections  at  any  point,  in  order  to  satisfy  the 
resonance  condition  the  impedance  looking  in  one  direction  should  be  equal  to  the 
negative  of  the  impedance  looking  in  the  other  direction.  In  our  analysis  first  we 
assume  that  the  two  ends  of  the  periodic  structure  are  open  circuited  and  since  we 
transform  all  the  impedances  to  the  one  end,  the  device  should  oscillate  at  the  frequency 
where  the  imag’  :ary  part  of  the  input  admittance  Yjn  becomes  zero.  The  conductance 
of  the  load  is  the  negative  of  the  real  pan  of  Yj^  (Imaginary  ofY =  0.0)  Therefore  by 
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connecting  a  load  with  resistance  to  one  end  of  the  parallel  plate  waveguide  the 

^  in 

oscillaton  condition  is  satisfied. 

12  Analysis  and  Results 

As  an  illustration,  a  parallel-plate  waveguide  loaded  with  four  QWITT  diodes 
has  been  chosen.  The  dimensions  of  the  waveguide  and  the  QWITT  diode  parameters 
designed  to  operate  at  94  GHz  are  as  follows: 

Width  of  the  waveguide  (a)  =  254  (im 
Thickness  of  the  waveguide  (b)  =  25.4  |im 
Distance  between  the  diodes  (d)  =  866  P-m 
Diameter  of  the  metal  post  and  the  diode  =  25.4  pm 
The  waveguide  is  assumed  to  be  filled  with  polyimide  (£^=3.4)  and  the  waveguide 
terminations  are  assumed  to  be  open  circuited  [6].  The  physical  parameters  of  the 
waveguide  should  permit  monolithic  circuit  fabrication  technology.  The  calculated 
input  admittance  as  a  function  of  frequency  is  shown  in  Fig.(2-7). 


Fig.(2-7):  Imaginary  part  of  the  input  admittance  versus  frequency 
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Fig.(2-8):  Dependence  of  the  oscillation  frequency  to  the  structure 
periodicity 

Maintaining  all  diode  and  circuit  parameters  the  same  as  those  at  94  GHz,  the 
distance,  d,  between  the  diodes  is  varied  to  obtain  other  resonant  frequencies.  At  any 
desired  operating  frequency,  the  distance  between  the  diodes  must  be  equal  to  half  the 
guide  wavelength.  Thus  as  the  oscillation  frequency  is  increased,  the  the  distance 
between  the  diodes  deceases  as  Vf  (Fig.2-8).  The  analysis  also  shows  how  variation  in 
the  width  and  thickness  of  the  waveguide  affects  the  oscillation  frequency  of  the  circuit 
designed  to  operate  at  94  GHz.  A  deviation  in  width  or  thickness  of  the  waveguide  by 
even  a  factor  of  four  causes  less  than  1%  change  in  the  oscillation  frequency.  Hence, 
the  width  and  thickness  of  the  waveguide,  which  determine  its  characteristic 
impedance,  do  not  affect  the  c''cillation  frequency  significantly.  This  is  a  distinct 
advantage  from  a  circuit  fabrication  point  of  view. 


Chapter3 

LARGE  SIGNAL  ANALYSIS 


The  optimum  design  of  microwave  circuits  containing  solid  state  devices 
requires  an  accurate  technique  for  predicting  their  nonlinear  performance.  Designs 
based  on  small-signal  parameters  can  roughly  estimate  the  frequency,  but  not  the 
power,  the  harmonic  content,  or  the  efficiency.  There  are  many  methods  of  nonlinear 
analysis,  each  one  having  its  own  specific  features  which  should  be  considered  based 
on  the  circuit  to  be  analyzed  and  the  desired  results.  Some  methods  are  better  for 
transient  analysis,  others  for  steady  state  response,  some  give  faster  results  and  some 
give  more  exact  results  but  are  very  slow. 

3.1  -  Different  Large  Signal  Analysis  Schemes: 

The  analysis  methods  may  be  divided  into  four  categories  [7] 

1)  Time  domain  and  direct  integration  methods:  for  transient  and  steady  state 
analysis. 

2)  Hybrid  analysis  using  time  and  frequency  domain  analysis  together:  for 
steady  state  analysis. 

3)  Nonlinear  transfer  functions:  for  weakly  nonlinear  circuit  analysis. 

4)  Voltage  update  method:  for  steady  state  analysis  of  a  nonlinear  network 
excited  by  a  sinusoidal  source. 

3.1.1  -  Time  Domain  Analysis: 

Time  domain  analysis  requires  solving  a  set  of  nonlinear  integro-differential 
equations.  Numerical  integration  techniques  can  be  used  to  solve  a  nonlinear  system  of 
equations  representing  a  dynamical  system  by  applying  appropriate  initial  conditions. 
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If  the  steady  state  solution  is  desired,  numerical  integration  should  be  performed  over 
many  cycles  till  the  transients  die  out.  If  the  solution  contains  slow  transients,  the 
computation  time  can  be  overwhelmingly  long.  This  is  often  the  case  for  high  Q 
circuits  such  as  cavity  stabilized  oscillators.  In  general  the  time  domain  approach 
requires  extensive  computation  time. 

3.12  -  Hybrid  Analysis:  Harmonic  Balance  Method 

In  the  harmonic  balance  method  [8]  the  periodic  or  quasi  periodic  steady  state 
solution  is  obtained  by  solving  the  nonlinear  differential  equations  representing  the 
circuit  in  the  frequency  domain.  This  is  carried  out  by  considering  as  unknowns  of  the 
analysis  problem  the  harmonic  components  Vj^,  1,^  of  the  Fourier  series  which  are  used 
to  expand  all  the  variables  in  the  circuit  (e.g.  the  voltages  v  and  the  currents  i)  Eq.(3-1). 

m 

v(i)  =  £  ReEVtA'] 

k^O 


m 


i(t)  =  Z 

k=0 


Re[lke^“k‘] 


(3-1) 


where  C)\  is  the  angular  frequency  of  the  spectral  component  and  m  is  the  number 
of  harmonics  needed  for  an  adequate  approximation  of  the  solution.  An  optimization 
algorithm  is  used  to  adjust  the  coefficients  of  the  Fourier  series  such  that  the  system  of 
equations  is  satisfied  with  least  square  error.  Although  the  harmonic  balance  method 
avoids  the  computationally  expensive  pioress  jf  numerically  integrating  the  dynamic 
equations,  its  main  disadvantage  is  the  large  number  of  variables  that  must  be 
optimized. 
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3.U  -  Volterra  Series  Approach 

When  the  circuit  to  be  analyzed  is  weakly  nonlinear,  an  approximated  steady- 
state  solution  can  be  obtained  by  Volterra  series.  The  input-output  relation  for  such 
nonlinear  circuits  can  be  viewed  as  a  generalization  of  the  convolution  integral  theorem 
in  the  analysis  of  linear  circuits.  If  x(t)  is  the  input  to  a  linear  system,  h(t)  is  its  impulse 
response  and  y(t)  is  its  output,  then: 

t  t 

Y(t)  =  J  X(t)  h(t-x)dx  =  J  X(t-x)h(x)dx  (3-2) 


In  the  case  of  a  nonlinear  system  the  input-output  behavior  can  be  described  by 
the  Volterra  functional  expansion  [9]. 

Y(t)  =  ho(t)  +  J  hi(t-Xi)X(Xi)dXi  ...+  1)1  h„(t-Xi,t-X2,...t-X„) 

X(Xi)X(X2)...X(x„)  dXjdX2....dXn  (3-3) 

where  hn(t,ti,....tn)  is  called  the  nonlinear  impulse  response  of  order  n.  Its  Fourier 
transform  Hn(fi,f2,...,fn)  can  then  be  called  the  nonlinear  transfer  function  of  order  n. 
In  order  to  analyze  a  nonlinear  circuit  by  means  of  Volterra  series,  it  is  necessary  to 
calculate  the  nonlinear  transfer  functions  Hi,H2,...,Hn.  As  in  any  series  expansion, 
this  approximation  is  useful  when  the  number  of  terms  required  for  an  adequate 
approximation  to  y(t)  is  not  very  large.  This  is  why  the  Volterra  series  approach  is 
often  used  to  analyze  the  systems  with  a  weak  nonlinearity. 
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3.1.4  -  Voltage  Update  Method 

This  method  is  applicable  to  the  analysis  of  a  circuit  containing  linear  and 

nonlinear  elements  which  is  excited  by  a  sinusoidal  source  E(t)  [10].  Like  the  harmonic 

balance  method  the  circuit  is  devided  into  two  parts;  one  contains  nonlinear  elements  (as 

well  as  some  linear  elements  if  necessary)  and  the  other  part;  contains  the  remaining 

linear  elements  and  the  exciting  source.  The  analysis  starts  with  selection  of  an 

estimated  value  of  voltage  across  the  terminal  connecting  the  linear  and  the  nonlinear 
section  VQ(t);  this  value  is  usually  taken  at  the  exciting  frequency,  without  harmonics. 

N  N  1  N 

Using  VQ(t),  the  current  in  the  nonlinear  circuit  IqCO  is  found.  Putting  I^Ct)  =  -  Io(t). 

the  fast  Fourier  transform  can  be  used  to  obtain  l|)(co). which  is  applied  to  the  linear 
circuit  VqCco)  is  found  by  superposition.  Then  using  an  inverse  Fourier  transform  we 
form  VqCO  and  compare  it  with  initial  Vq  (t).  If  the  difference  is  significant,  the 

iteration  proceeds  with  a  new  estimate  Vj(t),  and  continues  until  the  difference  is 

significant. 

3.2  -  The  Piecewise  Harmonic  Balance  Method  For  Large  Signal  Analysis 

Although  the  harmonic  balance  method  avoids  the  computationally  expensive 
process  of  numerically  integrating  the  dynamic  equations,  its  main  disadvantage  is  the 
large  number  of  variables  that  must  be  optimized.  If  the  system  of  equations  contains  n 
state  variables  and  each  of  them  requires  2m+l  Fourier  coefficients  (for  the  dc 
component  and  m  harmonics),  then  there  will  be  n(2m+l)  variables  to  be  adjusted  by 
the  optimization  routine.  This  makes  the  method  impractical  fc«-  large  systems,  unless  a 
way  is  found  to  reduce  the  number  of  optimized  variables,  in  order  to  reduce  the 
complexity  of  the  problem,  a  piecewise  harmonic  balance  method  [2]  can  be  used.  This 
method  is  essentially  an  harmonic  balance  technique  with  the  advantage  that  it  takes  into 
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account  the  fact  that  considerable  part  of  the  network  is  usually  linear.  In  this  method 
not  all  the  electrical  variables  in  the  circuit  (e.g.  nodal  voltages  and/or  branch  currents), 
but  only  those  relative  to  the  nonlinear  components  are  dealt  with  as  unknowns.  This  is 
used  to  reduce  the  number  of  optimized  variables. 

The  piecewise  harmonic  balance  method  may  be  described  by  the  following 

steps; 

1)  The  nonlinear  network  under  consideration  is  decomposed  to  a  minimum 
possible  number  of  linear  sub-networks  and  a  minimum  number  of  nonlinear 
ones.(Fig.3-l). 

2)  The  terminals  of  the  linear  sub  networks  are  excited  by  periodic  sources  in 
the  form  of  Fourier  expansions. 

3)  Through  the  solution  of  the  linear  sub-networks  in  the  frequency  domain  and 
by  a  suitable  optimization  routine,  the  coefficients  of  the  Fourier  series  are 
adjusted  so  that  the  equations  representing  both  the  linear  and  the  nonlinear 
sub- networks  are  satisfied. 

After  partitioning  the  circuit  into  linear  and  nonlinear  sub-netwoiks  as  shown  in 
Fig.  (3-1)  (i.e.  by  putting  only  the  intrinsically  nonlinear  components  into  the  nonlinear 
blocks),  only  electrical  variables  at  the  connection  ports  have  to  be  considered  as 
unknowns  of  the  harmonic  balance  problem. 
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The  linear  sub-network  response  at  the  connection  ports  is  obtained  in  the 
frequency  domain  by  the  set  of  equations: 

AjjVjj  +  =  Q  for  k=0, . ,m  (3-4) 

where  Vj^  and  are  the  column  vectors  of  the  k*  spectral  components  of  the  n  voltages 
and  currents  at  the  connection  ports;  Aj^  and  are  determined  by  the  linear  network 
and  Cif  represents  the  excitation  terms. 

The  nonlinear  sub-network  can  be  described  by  n  nonlinear  time-domain 
equations.  In  the  case  of  a  simple  resistive  voltage-controlled  nonlinear  elements  it  can 
be  written  as: 


i=f(v) 


(3-5) 
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or,  in  the  case  of  a  nonlinear  element  for  which  charge  storage  effects  should  also  be 
taken  into  account,  it  takes  the  form 


i=f(v,v) 


(3-6) 


where  v’  is  the  time  derivative  of  voltage.  In  the  more  general  case  for  complex 
nonlinear  devices,  the  nonlinear  equation  has  the  form 


f(v,v,...,v,i,i,...,i)  =  0 


(3-7) 


The  nonlinear  analysis  consists  of  searching  for  those  Vj^,  Ij^  which  satisfy  (by 
taking  into  account  Fourier  series  Eq.  3-1)  both  the  linear  frequency  domain  Eq.(3-4) 
and  the  nonlinear  time  domain  Eq.(3-7).  Once  the  harmonic  balance  problem  has  been 
solved,  all  the  other  variables  of  the  circuit  can  be  determined  by  solving  linear  sub¬ 
network  equations  in  the  frequency  domain.  Eq.(3-4)  can  be  written  in  the  form 

Vk  =  Zklk  +  for  k=0, . m  (3-8) 


By  transforming  this  into  time  domain  and  substituting  it  into  Eq.(3-7)  we  get  a  set  of 
time-dependent  equations 


gvlo» . Ik’ 


,i„,o-o  for-y^sts-y^ 


(3-9) 
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where  T  is  the  period  of  the  circuit  response.  In  the  case  of  quasi-periodic  operation  T 
goes  to  infinity. 

Since  the  variables  v(t),  i(t)  have  been  approximated  by  a  Fourier  series  with  a 
finite  number  of  terms  (m),  it  is  not  generally  possible  to  find  a  set  of  harmonic 
components  Iq,. of  currents  which  exactly  satisfy  Eq.(3-9)  for  every  t,  but 
only  an  approximate  solution  which  gives  the  minimum  discrepancy  in  Eq.(3-9).  The 
minimum  discrepancy  criterion  can  be  defined  in  different  ways,  and  in  particular ,  by 
the  least-square  error  function. 

n  T/2 

EFdo, . ,Ik, . I^)  =  £  T  Jjg  n(Io . Ik . Im.t)  dt  (3-10) 

Therefore  the  best  approximated  solution  to  the  harmonic  balance  problem  can  be 
computed  by  finding  harmonic  components  Iq,.. of  the  currents  which  minimize 
EF  by  means  of  a  numerical  optimization  algorithm.  It  should  be  noted  that  we  could 
optimize  coefficients  of  Fourier  series  for  voltage  wave  forms  at  each  interconnection 
instead  of  current  wave  forms  as  in  the  above  case. 

In  the  analysis  of  oscillators  ,  it  is  also  necessary  to  consider  as  unknown  the 
period  of  the  circuit  response.  In  this  case  the  unknown  variables  of  the  optimization 
problem  are  not  only  lo,....,!^,  but  also  the  period  T. 

When  the  least  square  approximation  criterion  is  used  to  define  the  error 
function  ,  a  general  purpose  numerical  optimization  algorithm  such  as  a  quasi  Newton 
gradient  method  can  be  used  to  minimize  the  error  function. 

In  the  harmonic  balance  technique  there  is  always  a  need  for  a  suitable  initial 
estimate  for  the  solution.  This  initial  estimate  strongly  affects  the  convergence 
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characteristics  of  the  optimization  procedure.  Since  in  practice  it  is  very  difficult  to 
estimate  the  unknown  harmonic  components  of  the  variables,  convergence  problems 
would  be  encountered  in  strongly  nonlinear  circuits,  unless  suitable  strategies  were 
adopted.  One  way  is  to  start  from  small  values  of  the  excitation  terms,  so  that  nonlinear 
effects  are  not  strong,  then  increase  their  amplitude  gradually  in  successive  steps,  till 
the  nominal  values  are  reached.  At  every  step  the  solutions  found  in  previous  steps  are 
used  as  initial  estimate  for  the  iterative  analysis  algorithm.  (It  should  be  considered  that 
the  above  can  not  be  applied  to  the  case  of  oscillator  analysis.) 

Since  we  are  interested  only  in  the  steady  state  analysis  of  microwave 
oscillators,  the  harmonic  balance  seems  to  be  the  proper  method  of  analysis.  For  the 
optimization  part  a  general  purpose  optimization  routine  from  the  IMSL  library 
(UMINF)  can  be  used.  The  subroutine  NONLINEAR  in  our  program  contains  the 
nonlinear  device  model.  By  modifying  this  subroutine  different  devices  can  be  used  as 
the  nonlinear  part  of  the  oscillator  circuit.  In  order  to  test  the  convergence  behavior  of 
the  nonlinear  analysis  program  two  different  test  problems  are  solved.  This  helps  us  to 
investigate  the  performance  of  the  optimization  routine  and  study  the  sensitivity  of 
convergence  to  the  initial  guess. 

32.1  -  Test  I:  Van  der  Pol  Oscillator 

Van  der  Pol's  equation  seems  to  be  a  classical  example  of  a  nonlinear  oscillatory 
system.  It  is  of  the  form 

Xi  =  X2 

X2  =  p(l-Xj)X2-Xi 


(3-11) 
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The  network  realization  of  Van  der  Pol's  equation  is  shown  in  Fig.(3-2) 


Fig.  (3-2):  Circuit  representation  of  the  Van  der  Pol's  oscillator 

The  factor  |i  determines  the  amount  of  nonlinearity.  If  |j,  is  very  small,  the  term  p.(l  - 
Xj)  can  be  neglected.  In  this  case  the  differential  equation  represents  a  simple 

harmonic  oscillatory  system  which  consists  of  a  1  F  capacitor  in  parallel  with  a  1  H 
inductor.  In  this  case  the  period  of  oscillation  is  2n.  If  on  the  other  hand  |j,  is  not 
small,  because  of  nonlinearity  effects,  the  oscillation  waveform  will  contain  higher 
harmonics.  Harmonic  contents  cause  the  period  of  oscillation  to  increase  slightly.  This 
is  because  the  voltage  is  not  a  single  harmonic  any  more.  Acting  with  the  nonlinear 
element,  it  produces  a  current  which  is  far  from  a  simple  sinusoid.  One  component  of 
this  current  is  at  the  fundamental  frequency.  The  frequency  must  be  such  that  the 
fundamental  current  multiplied  by  the  impedance  of  the  LC  circuit  at  this  frequency 
gives  the  fundamental  voltage.  This  condition  requires  that  the  resonance  frequency  of 
the  oscillation  be  slightly  lower  than  the  resonance  frequency  of  the  LC  circuit.[l  1]  In 
this  analysis  we  choose  |i=l  (strong  nonlinearity). 

In  order  to  analyze  the  circuit  using  the  harmonic  balance  method,  the  circuit  is 
divided  into  nonlinear  and  linear  subcircuits.  The  nonlinear  subcircuit  is  shown  by  the 
dashed  lines  (Fig.3-2).  For  the  analysis  we  assume  the  linear  and  the  nonlinear 
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subcircuits  are  excited  by  a  voltage  source  V(fj„.fj„)  which  contains  m  harmonics.  For 
the  steady  state  resonance  condition  to  be  satisfied,  the  impedance  looking  into  the 
nonlinear  element  should  be  equal  to  the  negative  of  the  impedance  seen  from  the  linear 
subcircuit  port.  Therefore  if  the  voltage  V  is  equal  to  the  true  steady  state  voltage  which 
should  appear  across  the  linear  and  nonlinear  subcircuit's  terminals,  the  sum  of  the 
currents  flowing  into  the  linear  and  nonlinear  subcircuits  should  be  equal  to  zero. 

The  nonlinear  analysis  program  requires  an  initial  guess  (Vg)  for  the  voltage  V. 
Generally  the  initial  guess  contains  an  approximate  value  for  the  frequency  and 
amplitude  of  the  fundamental  oscillation  (it  can  also  contain  amplitude  and  phase  of  the 
higher  harmonics).  Then  the  current  flowing  into  the  linear  circuit  (Iimear(fi.-fm').  ni' 
can  vary  between  1  to  m  depending  on  the  number  of  harmonics  in  the  initial  guess)  is 
calculated  in  the  frequency  domain.  After  sampling  the  voltage  Vg  in  the  time  domain, 
the  current  through  the  nonlinear  element  is  determined  (Inonimear(t)  )•  Through  an 
iterative  procedure  using  the  optimization  routine  frequency,  amplitudes  and  phases  of 
the  n  harmonics  of  the  voltage  (Vg)  can  be  chosen  such  that  an  error  function  defined  as 
r.m.s  value  of  the  sum  of  currents  through  the  linear  and  the  nonlinear  subcircuits  is 
minimized.  The  voltage  waveform  Vg  that  minimizes  the  error  function  is  the  steady 
satate  solution  to  the  problem.  Since  the  current  through  the  linear  circuit  is  in  the 
frequency  domain  and  the  current  through  the  nonlinear  circuit  is  in  the  time  domain, 
one  can  define  the  error  function  in  the  time  (by  transferring  Iiinear(f!"fm')  to  the  time 
domain)  or  in  the  frequency  domain  (by  transferring  InoniinearCt)  to  the  frequency 
domain  using  an  FFT  algorithm).  The  error  functi**  n  in  the  time  domain  is  given  by 

EF,=  I  ( W„(t)  +  )  2  (3-12) 
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where  T  is  the  period  of  oscillation  and  t|,t2.— are  the  sampling  points. 

The  error  function  in  the  frequency  domain  is  given  by 

m 

EFf  =  s  (  Wf,)  +  ^  (3-13) 

1=1.2.. 

where  m  is  the  number  of  harmonics  considered  in  the  analysis.  Either  of  the  two  error 
functions  can  be  used  in  the  optimization  procedure,  but  as  presented  below,  different 
performance  is  observed  with  the  two  different  error  functions. 

The  sensitivity  to  the  initial  guess  of  the  program  using  Erfj  is  less  even  though 
the  program  takes  a  longer  time  to  converge.  Results  of  several  runs  containing  the 
initial  guess  and  the  final  values  after  convergence  are  given  below.  After  each  result 
the  wave  form  across  the  nonlinear  device’s  terminals  is  plotted.  The  period  and 
amplitude  of  the  waveform  is  very  close  to  the  one  calculated  by  Nakhla  [2]. 

Initial  guess: 

amplitude  of  the  first  harmonic:  2.0 
phase:  0 

period  of  oscillation:  5.0 

Final  value: 

period  of  oscillation:  6.6645s 
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The  nonlinear  analysis  program  which  uses  the  error  function  in  the  frequency 
domain  is  less  sensetive  to  the  initial  guess. 

-  Test  2:  Diode  With  Negative  Conductance 
To  test  the  nonlinear  analysis  program  again  the  large  signal  analysis  of  a 
negative  resistance  diode  oscillator  is  performed.  The  negative  resistance  diode  is 
inserted  in  a  parallel  plate  waveguide  via  an  inductive  post 

The  diode  has  a  negative  dynamic  conductance  which  is  obtained  by  application 
of  a  proper  bias  voltage.  Some  examples  of  negative  conductance  diodes  are  IMP  ATT, 
Tunnel,  TRAPATT  and  QWITT  diodes.  As  a  first  approximation,  the  frequency 
dependence  of  the  negative  conductance  may  be  neglected  as  long  as  the  frequency 
range  over  which  the  diode  operates  is  small  enough.  The  negative  conductance  may 
be  represented  in  the  frequency  range  which  the  above  approximation  is  valid,  by  a 
function  relating  high  frequency  current  to  voltage.  This  function  is  nonlinear  for  large 
signals  applied  across  the  diode's  terminals.  This  nonlinear  relation  may  be 
approximated  by  a  cubic  function  shown  in  Fig.(3-5).[12] 

If  the  operation  point  is  assumed  to  be  the  origin  of  the  coordinate  system,  the 
instantaneous  value  of  the  current  given  by: 


I  =  -ajV  +  a3V3 


(3-14) 
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Fig.(3-5):  cubic  approximating  curve 

The  dynamic  conductance  is  the  derivative  of  the  current  with  respect  to  the  voltage.  It 
is  given  by: 

G  =  -ai  +  3a3V2  (3.15) 

This  shows  that  as  the  instantaneous  value  of  the  voltage  across  the  diode 
increases  the  dynamic  negative  conductance  of  the  diode  decreases.  The  equivalent 
circuit  of  the  negative  conductance  diode  is  shown  in  Fig.  (3-6) 


Fig.(3-6):  Equivalent  circuit  of  a  negative  resistance  diode 


The  negative  conductance  (Gj)  is  shunted  by  the  junction  capacitance  (0^).  and 
are  series  inductance  and  resistance  respectively. 
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Both  the  diode's  negative  conductance  and  junction  capacitance  are  assumed  to 
be  frequency  independent.  The  diode's  negative  conductance  may  be  expressed  as 


G<,=  -G„(l-av5 


(3-16) 


where  Gd  is  the  r.m.s  negative  conductance  of  the  diode,  Gq  is  the  negative  value  of 
the  diode  small  signal  conductance,  a  is  a  diode  parameter  which  can  be  determined  by 
measurements  on  an  oscillating  diode  and  V©  is  the  r.m.s  value  of  the  voltage  across 
the  diode.  The  comparison  of  Eq.(3-14)  and  Eq.(3-16)  results  in 

Go  =  ai  (3-17) 

a  =  (3-18) 


The  level  dependence  of  the  junction  capacitance  may  be  approximated  by  an 
expression  similar  to  that  giving  the  level  dependence  of  the  conductance. 

Cd  =  Cdo(l+PVo2)  (3-19) 

Where  P  is  a  diode  parame  er.  The  impedance  across  the  diode's  terminals  is  given  by: 
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G^  +  jmq  (3-20) 

<35  + 

The  impedance  is  dependent  on  the  amplitude  of  the  voltage  across  the  diode  because 
Gd  and  Cd  are  dependent  on  it.  Due  to  the  junction  capacitance  and  series  resistance, 
the  negative  conductance  of  the  diode  decreases  with  inceasing  frequency  even  though 
the  junction  conductance  Gd  is  assumed  to  be  independent  of  frequency. 

322  Ji  -  Estimation  of  the  Amplitude  and  Frequency  of  Oscillation 
To  provide  resonable  estimate  of  the  frequency  and  the  amplitude  of  oscillation 
to  be  used  as  initial  guess  for  the  harmonic  balance  method  the  following  simplified 
approach  is  performed. 

The  simplified  equivalent  version  of  a  negative  conductance  oscillator  is  shown 
in  Fig.  (3-7) 


X 


Fig(3-7):  A  negative  resistance  oscillator 


This  study  is  based  on  the  assumption  that  only  the  fundamental  frequency 
appears  across  the  diode's  terminals  and  all  the  higher  harmonics  are  shorted  out. 
Through  this  analysis  we  can  determine  an  approximate  value  for  the  oscillation 


firequency  and  the  amplitude  of  oscillation.  We  shall  use  this  approximate  solutions  as 
the  initial  guess  for  the  harmonic  balance  analysis. 

In  the  stable  state,  the  admittance  across  the  terminals  x-x  in  Fig.(3-7)  seen  in 
one  direction  should  be  equal  to  the  negative  of  the  admittance  seen  in  the  other 
direction. 

The  nonlinear  admittance  of  the  diode  is  given  by 


Yd=  -God  -  avj  +  jtoCdod  + 

(3-21) 

where  Vj  is  r.m.s  value  of  the  voltage  across  the  diode's  terminals. 

requirements  for  steady  state  oscillation: 

To  satisfy  the 

Oo(l-av5)  =  G, 

(3-22) 

-OjQiod  +  PV^  =  (oCp  -  1/  (oLp 

(3-23) 

The  power  generated  by  the  diode  is  given  by 

P  =  Giv2  =  Goa-av2)v2 

(3-24) 

In  order  to  get  the  maximum  output  power,  we  differentiate  the  Eq.(3-24)  with  respect 

2 

to  Vj  and  equate  it  to  zero  to  get 
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a  = 


(3-25) 


where  Vm  is  the  rms  voltage  that  yields  the  maximum  power.  Substituting  Eq.(3  25) 
into  Eq.(3-24)  we  get  an  expression  for  maximum  output  power. 


Pm 


(3-26) 


The  normalized  power  can  be  expressed  as 


GiV 


Pn  =  ^= 

“m 


d  „  2 

=  2giVd 


(3-27) 


where 


Gl 

gl  =  ^ 


and 


vd  =  Y^ 


m 


Using  Eqs.  (3-22)  and  (3-25)  we  get 


1  ^  2 
gl=l  - 


(3-28) 

(3-29) 


(3-30) 


By  combining  Eqs.  (3-27)  and  (3-30)  we  get 
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Pn=4gi(l  -  gi) 


(3-31) 


This  equation  shows  the  dependence  of  the  power  on  tne  load.  The  highest  power  is 
obtained  when  gi  =^or 


(3-32) 


This  shows  in  order  to  get  the  maximum  power  the  conductance  of  the  load  should  be 
half  the  small  signal  conductance  of  the  diode. 

The  oscillation  frequency  may  be  obtained  from  Eq.(3-23).  If  the  amplitude  of 
voltage  across  the  diode's  terminals  is  small. 


^Lp(Cdo  +  Cp) 


(3-33) 


The  oscillation  frequency  is  decreased  with  increasing  voltage  across  the  diode. 

3  J  J.b  -  The  Harmonic  Balance  Results 

A  sample  problem  is  designed  to  test  the  nonlinear  analysis  program.  A 
negative  conductance  diode  is  inserted  in  a  parallel  plate  waveguide.  The  I-V 
characteristic  for  the  test  diode  is 


I  =  O.OKV  -  3.0)3  +  8.56x10-2  (V  -  3.0)  -h3.0 


(3-34) 


where  ai=8.56xl0'2  and  a3=0.01, therefore  the  small  signal  negative  conductance  Go  = 
8.56x10'2.  In  the  test  case  it  is  assumed  that  the  junction  capacitance  Cd  is  frequency 
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and  level  independent.  The  value  of  Cd  is  chosen  to  be  equal  to  0.1  pF.  For  the 
analysis  the  Cd  is  embedded  in  the  linear  part  of  the  circuit.  By  comparing  Eq.(3-24) 
with  Eqs.(3-17)  and  (3-18)  a  =  0.175  .  In  order  to  obtain  maximum  power  the 
resistance  of  the  load  seen  by  the  diode  should  be  equal  to  23.4  Q.  Based  on  the  small 
signal  analysis  the  following  oscillator  circuit  is  designed  to  oscillate  at  frequency  75.3 
GHz. 


negative  resistance  device 


Fig(3-8):  A  negative  resistance  oscillator 


b(height)=.o254 

a(width)=.254 

d=.815 


Zj  =20.45  ohms 
Ik;, c= 10.0  ohms 


V 


bias 


=33.0  volts 
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Based  on  the  simplified  nonlinear  analysis  that  was  explained  in  the  last  section 
the  initial  guess  for  the  amplitude  of  oscillation  in  the  above  case  is  calculated  to  be 
about  2.4  volts  peak,  (using  Eq.  3-25)  Also  the  resonance  frequency  of  the  linear 
circuit  is  75.3  GHz. 

After  performing  the  harmonic  balance  analysis  on  the  above  circuit  the  foilwing 
results  are  obtained. 

Number  of  harmonics  considered  across  the  diode's  terminals:  5 
Initial  guess  for  the  oscillation  frequency:  75.3  GHz 
Initial  guess  for  the  amplitude  of  the  fundamental  component:  2.4  volts 
Results  after  convergence: 

Vdc  =  3.0  volts 

Frequency  of  oscillation  =  74.39  GHz 
Voltage  components  across  the  diode's  terminals: 


N 

1 

2 

3 

4 

5 


Amplitude  (volts) 

2.4137 

2.99x10^ 

0.1746 

7.649x10^ 

3.72x10'^ 


Phase(radians) 

0.0 

1.74 

1.98 

-2.07 

-2.48 


It  can  be  seen  that  due  to  generation  of  higher  harmonics  ‘ne  oscillation 
frequency  is  lower  that  the  resonance  frequency  of  the  linear  circuit.  The  frequency 
shift  is  small  (about  1.2%).due  to  the  fact  that  it  is  a  second  order  effect.  Also  the 
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amplitude  of  the  fundamental  frequency  of  oscillation  is  very  close  to  the  expected 
value. 

To  check  the  convergence  of  the  analysis  we  chose  a  different  initial  guess  and 
mn  the  program  again. 

Number  of  harmonics  considered  across  the  diode's  terminals:  5 
Initial  guess  for  the  oscillation  frequency:  76.0  GHz 
Initial  guess  for  the  amplitude  of  the  fundamental  component:  6.0  volts 
Results  after  convergence: 

Vdc  =  3.003  volts 

Frequency  of  oscillation  =  74.395  GHz 
Voltage  components  across  the  diode's  terminals: 


N 

1 

2 

3 

4 

5 


Amplitude  (volts) 

2.4132 

9.0x10^ 

0.17487 

2.659x10’^ 

4.06x10'^ 


Phase(radians) 

0.0 

-0.511 

2.088 

-2.185 

-2.324 


The  initial  guess  for  the  amplitude  of  oscillation  is  bigger  than  twice  the  final  value. 
Convergence  to  the  same  amplitude  occurs  with  the  initial  values  chosen  to  be  smaller 
than  the  tn>e  value.  This  proves  the  existance  of  a  stable  limit  cycle  with  an  amplitude 
of  2.41  volts.  Of  course  if  the  value  of  the  initial  guess  is  too  small  then  the  amplitude 
will  converge  to  zero  which  is  also  a  stable  focus  (No  oscillation).  Unfortunately  the 


harmonic  balance  method  does  not  give  us  any  insight  about  the  build  up  of  oscillation, 
because  that  is  a  transient  phenomenon.  In  the  design  of  oscillators  one  should  be  sure 
that  the  oscillation's  start  up  conditions  are  satisfied  [12]. 

33  -  Large  Signal  Analysis  of  A  QWTTT  Oscillator 

Using  the  harmonic  balance  method  the  large  signal  analysis  of  the  QWITT 
diode  oscillator  is  performed.  The  analysis  is  similar  to  the  one  in  test  #2.  The  only 
difference  is  that  the  subroutine  NONLINEAR  contains  the  large  signal  model  of  the 
QWITT  diode.  As  before  the  QWITT  diode  is  inserted  in  a  parallel  plate  waveguide  via 
an  inductive  post.  The  nonlinear  analysis  is  performed  considering  only  one  harmonic 
across  the  QWITT  diode.  Due  to  the  fact  that  the  amplitudes  of  the  higher  harmonics 
are  much  smaller  than  the  amplitude  of  the  fundamental  frequency  of  operation, 
neglecting  higher  harmonics  will  not  have  much  effect  on  the  accuracy  of  our  analysis. 
The  large  signal  analysis  is  performed  only  for  a  single  diode  oscillator. 

33.1  -  Large  Signal  Model  of  The  QWITT 

We  want  to  find  the  current  through  the  QWTTT  diode  as  a  function  of  dc  bias, 
ac  voltage  at  the  fundamental  frequency  across  the  diode,  and  frequency  of 
operation.. The  large  signal  diode  model  [13]  uses  a  piecewise  linearized  fit  to  an 
experimental  dc  IV  characteristic  of  the  quantum  well  to  describe  current  injection. 
Using  the  Ramo- Shockley  theorem  for  transit  time  effects,  a  numerical  integration  of 
the  injected  charge  pulse  over  the  entire  rf  cycle  is  performed  to  obtain  the  device 
current.  The  electric  field  and  the  carrier  velocity  in  the  drift  region  are  assumed  to  be 
constant.  Static  quantum  well  I  -  V  characteristics  are  assumed  to  be  valid  at  all 
frequencies.  Fourier  analysis  of  the  device  current  gives  the  device  conductance  and 
susceptance  per  unit  area. 
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First  we  need  to  find  an  expression  for  the  total  current  density  as  a  function  of 
injection  angle  0.  Consider  that  a  pardcle  is  injected  at  an  angle  0  from  the  quantum 
well  device  into  the  drift  region.  If  we  assume  that  the  electric  field  in  the  drift  region  is 
strong  enough  (high  bias  voltage)  then  the  particle's  velocity  in  the  drift  region  can  be 
assumed  to  be  equal  to  the  saturation  velocity  of  electrons  v^  in  that  region.  By  the  time 

the  particle  reaches  the  contact  at  the  end  of  the  drift  region,  its  phase  angle  will  be 
0+0^  where  0^  is  the  drift  angle  of  the  drift  region.  At  angle  0  all  particles  injected 

from  0-0^  up  until  0  are  still  drifting  through  the  device.  We  can  write 

Q(t)  =  f  ‘jQwdt’  (3-35) 

where 

Q(t)  =  Total  charge  drifting  through  the  device  at  time  t. 
jQWi(t)  =  Particle  ciurent  density  through  the  quantum  well  at  time  L 
It  is  assumed  that  Jq^=  JqwCVqw)  ^  known  J-V  characteristic  and  that  the  transient 

J-V  characteristic  is  approximately  equal  to  the  steady  state  J-V  characteristic. 

We  can  change  the  integration  variable 

2n  T 
so 

dt  =  — d0 
27rf 


where  f  is  the  frequency  of  operation.  Eq.(3-35)  can  be  written 
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Q(e)  = 


By  the  Ramo  -  Shockley  theorem,  we  can  write 


J(e)  = 


Q(e).v, 

w 


where 

Vj  =  saturation  velocity  of  electrons  in  the  drift  region 
W  =  drift  region's  length 
From  Eq,(3-38)  and  Eq.(3-39)  we  can  write 

J(0)=-^^  f  jQw(0‘)d0‘ 

27cfW  J  ^ 


27C  V 

and  since  —  W  =  and  X  =  -^  we  get 

X  ^ 


W  = 


0d^s 

27Cf 


therefore  by  substituting  Eq.(3-41)  into  Eq.(3-40)  the  following  is  obtained 


1(0)  =  ^ 

0d 


(3-38) 


(3-39) 


(3-40) 


(3-41) 


(3-42) 
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Given  Jqvv(Vqv^)  one  can  calculate  J(0)  by  numerical  integration  of 
Jqw(V qw(1  +  Ysin(0)))  where  it  is  assumed  that 


V  =  Vo(l  +7  sin(0))  =  total  voltage  across  device 


(3-43) 


QC 

Vqw=  ^qw(1  Y  an(0))  =  voltage  across  quantum  well 


(3-44) 


where 


Vq  =  the  dc  voltage  across  the  QWITT 
7  =  ac  modulation  factor 


Vi  =  Vo7sin(0) 


(3-45) 


where  Vj  is  the  ac  voltage  component  at  the  fundamental  frequency. 


Ji(Q)  =  hs  sin(0)  -I-  JicCos(0) 


(3-46) 


where  J^q  is  the  ac  component  of  current  at  the  fundamental  frequency.  We  can  find 
Jjs  and  Jjc  from  Fourier  analysis. 


Ju  =  -  sin(0)  d0 

7t«/o 


(3-47) 


Jic=-  r"‘j(0)cos(0)d0 

7C  •'o 


(3-48) 
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from  Eq.(3-42),  (3-47)  and  (3-48)  we  get 


sin(0)( 


2* 


cos(0)( 


(3-49) 


(3-50) 


33.2  -  Results  of  Harmonic  Balance  Analysis 

The  large  signal  conductance  and  susceptance  as  a  function  of  frequency  for  a 
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QWITT  diode  (10  |im  diameter)  operating  in  the  -  j-  mode  is  shown  in  Fig.(3-9).  As 

one  can  see  the  negative  conductance  of  the  diode  decreases  with  increasing  frequency. 
Also  the  susceptance  of  the  QWITT  which  is  due  to  the  parallel  capacitance  of  the  diode 
increases  as  frequency  increases.  These  two  effects  both  reduce  the  efficiency  of  the 
diode  at  very  high  frequencies,  [14] 

The  plot  of  the  output  power  of  a  QWITT  oscillator  as  a  function  of  frequency 
is  shown  on  Fig. (3- 10),  The  length  of  the  drift  region  is  optimized  to  generate 
maximum  power  at  each  frequency.  The  linear  part  of  the  circuit  was  designed  to 
match  the  diode  impedance  perfectly.  As  shown  the  power  decreases  with  frequency 


increase. 
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Frequency  (GHz) 


Fig.  (3-9);  Large  signal  conductance  and  susceptance  as  a  function  of 
frequency  for  a  QWITT  diode(10  mm  diameter)  operating  in  the-^  mode. 


10  lor  1000 

Frequency  (GHz) 

Fig.(3-10):  Output  power  as  a  function  of  frequency  for  a  single  QWITT  diode 


44 


0.0  0.1  0.2  0.3  0.4 

Gamma  O'ac/Vdc) 


Fig.(3-1 1):  Output  power  and  efficiency  as  a  function  of  gamma 


On  Fig.(3-1 1)  the  plot  of  output  power  and  efficiency  as  a  function  of  the  ratio 
of  the  ac  to  the  dc  voltage  (gamma)  across  the  diode  is  shown.  The  oscillation 
frequency  is  94  GHz.  At  each  gamma  point  an  optimum  load  is  used  for  maximum 
power  generation.  It  is  clear  that  there  is  an  optimum  ac  modulation  across  the  diode 
for  maximum  power  generation.  As  the  amplitude  of  the  ac  voltage  across  the  diode 
increases  the  impedance  of  the  diode  changes  (the  negative  conductance  decreases)  and 
that  causes  the  power  to  drop  from  its  optimum  value. 


Chapta*4 

EXPERIMENT 


Two  experiments  were  performed  to  demonstrate  the  oscillation  capability  of  the 
QWnr  at  microwave  frequencies.  For  both  experiments  QWTTT  diodes  with  different 
drift  region  length  were  used  to  determine  the  dependence  of  the  power  generated  on 
the  drift  region  lengths.  The  power  obtianed  at  different  frequencies  was  lower  than  the 
maximum  theoretically  predicted  power  for  the  QWITT  diode.  The  physical  structure 
of  the  device  is  shown  on  Fig.(4-1)  [15] 


0.5  urn  n  GaAs _ 

,  17  -3 

100  A  1x10  cm  GaAs 
16  .3 

100  A  1x10  cm  GaAs 
50  A  undoped  GaAs 
17  A  AlAs  barrier _ 

50  A  GaAs  well _ 

17  A  AlAs  barrier 
50  A  undoped  GaAs 

16  -3 

5x10  cm  doped  GaAs 

n  GaAs  buffer _ 

+ 

n  GaAs  substrate 


500A,  lOOOA,  and  2000A 


Fig.(4'l):  physical  structure  of  the  QWITT 


4.1  -  QWTTT  Diode  &i  A  Rectangular  Waveguide 

Three  different  QWTTT  diodes  with  drift  regions  500A,  lOOOA  and  2000 A  were 

7  2 

examined.  Devices  had  a  nominal  area  of  3.3-3.8x10'  cm  with  a  height  of  around 
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0.5  )xm.  The  devices  were  mounted  in  a  WG-16  (8.2  -  12.4  GHz)  waveguide  using  a 
micrometer  controlled  post  and  a  whisker  contact.  One  side  of  the  waveguide  was 
connected  to  a  tunable  short  (Fig.  4-2)  The  frequency  of  the  oscillator  could  be 
adjusted  by  varying  the  position  of  the  short  and  the  bias  voltage. 


Fig.(4-2):  waveguide  oscillator  using  QWTTT  diode 


The  rf  performance  of  three  diodes  is  given  in  table  (4-1).  The  output  powers  obtianed 

were  less  than  the  theoretically  predicted  values.  This  may  be  due  to  the  fact  that  the 

impedance  of  the  diode  was  not  matched  to  the  characteristic  wave  impedance  of  the  full 

height  waveguide.  The  maximum  power  was  obtained  using  a  device  with  the  drift 
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region  lengu  of  2000 A  which  is  very  close  to  the  mode  of  operation  at  10  GHz 
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Length  of  Drift 

e 

Region  (A) 

Specific  Negative  Resistance 
(10^  Ohm-cm^) 

Output  Power  (|iW) 

Oscillation  Frequency 
(GHz) 

500 

2.6 

3 

10.2 

1000 

4.9 

10 

10.4 

2000 

19 

30 

10.6 

Table  (4-1) 


42  -  Planar  Structure 

The  QWnr  diodes  were  mounted  on  a  planar  microstip  structure.  The 
maximum  power  obtained  was  about  900  |xW  at  4.2  GHz.  One  reason  for  obtaining 
higher  power  using  the  planar  circuit  is  that  the  impedance  of  the  diode  is  better 
matched  to  the  impedance  of  the  microstrip  structure.  As  in  the  case  of  the  waveguide 
oscillator,  here  the  maximum  pwwer  was  obtained  with  the  diode  with  2000A  drift 
region  length.  The  rf  performance  of  the  planar  structiue  is  given  in  table  (4-2). 


V 
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Length  of  Drift 

0 

Region  (A) 

Specific  Negative  Resistance 
(lO^Ohm-cm^) 

Output  Power  (M-W) 

Oscillation  Frequency 
(GHz) 

240 

10.0 

500 

2.6 

275 

7.0 

2000 

19 

910 

4.2 

Table  (4-2) 


> 


Chapta*5 

Conclusion 


We  have  applied  small  and  large  signal  models  for  the  QWl  Ti  diode  in  order  to 
prove  its  feasibility  for  use  in  millimeter  and  submillimeter  wave  oscillators.  Small 
signal  analysis  was  applied  to  a  periodic  power  combining  structure  to  determine  the 
design  parameters  and  oscillation  frequency  of  the  oscillatior.  The  large  signal  analysis 
for  a  single  QWITT  diode  oscillator  used  the  piecewise  harmonic  balance  method  to 
determine  the  precise  amplitude  and  frequency  of  oscillation.  The  calculations  showed 
that  the  QWITT  diode  is  capable  of  producing  about  1  mW  power  at  frequencies  up  to 
600  GHz.  Furthermore  an  experiment  was  performed  which  verified  the  capability  of 
the  QWITT  diode  for  the  generation  of  microwave  power. 

Future  experiments  are  planned  which  will  demonstrate  operation  at  higher 
power  and  higher  frequencies.  Also  a  periodic  power  combining  structure  will  be 
fabricated  to  increase  output  power  of  the  oscillator. 
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